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Abstract — In this paper, we consider the amplified-and-forward 
relaying in an OFDM system with unitary linear processing at 
the relay. We proposed a general analytical framework to find 
the unitary linear processing matrix that maximizes the system 
achievable rate. We show that the optimal processing matrix is 
a permutation matrix, which implies that a subcarrier pairing 
strategy is optimal. We further derived the optimal subcarrier 
pairing schemes for scenarios with and without the direct source- 
destination path for diversity. Simulation results are presented 
to demonstrate the achievable gain of optimal subcarrier pairing 
compared with non-optimal linear processing and non-pairing. 

L Introduction 

OFDM-based relaying combines the advantages of both 
OFDM and relaying techniques to improve network perfor- 
mance and efficiency. It is essential that the relaying tech- 
niques can be maximally explored in the OFDM system. In 
a narrow-band single carrier system, the relay retransmits the 
processed version of the received signal over the same carrier. 
For an OFDM system with multiple subcarriers, the relay has 
an additional frequency dimension available, which can be 
exploited to process the incoming signals adaptively based on 
subcarrier strength for relaying purpose, and thus can poten- 
tially improve the overall relay performance. Without such 
consideration, directly applying a relaying method optimized 
for single-carrier systems will not be optimal in an OFDM 
system. In this paper, we address how to optimally perform 
channel aware linear processing of the incoming signals over 
multiple subcarriers at the relay to maximize the relaying 
performance in an OFDM system. 

Subcarrier pairing, which maps incoming and outgoing 
subcarriers at the relay, can be viewed as a special case of 
unitary linear processing. It was first proposed independently 
in |1| and |2) for a single amplified-and-forward (AF) OFDM 
relay system. For relaying without the direct source-destination 
path available, 11] used integer programming to find the 
optimal pairing that maximizes the sum SNR, whereas from a 
system-design perspective, |2| proposed a subcarrier pairing 
scheme optimal in the noise-free case, assuming uniform 
power allocation. These works sparked interests for researchers 
in this area. Subsequently, a joint subcarrier pairing and power 
allocation scheme was studied in |3|, although the power 
constraint at the relay cannot be guaranteed in the presented 
scheme. Furthermore, the sorted subcarrier pairing was studied 



in MIMO OFDM systems H, as well as in decode-and- 
forward (DF) relay systems Q. All of the above works 
focus on the relay path only, without direct-path transmission, 
perhaps partially due to the difficulty of finding subcarrier 
pairing in such a case. For relaying with direct path available 
for full diversity, two suboptimal subcarrier pairing schemes 
were studied in |6|. 

In this paper, we propose a framework to address the optimal 
unitary linear processing for relaying in an OFDM system. A 
half-duplex two-hop AF relay network is considered, where 
a unitary matrix is applied to the incoming signals across 
subcarriers at the relay, and the relay then forwards the am- 
plified version of the processed signals. We propose a general 
framework that enables us to analyze the effect of unitary 
linear processing matrix on the achievable rate of such relaying 
system. Constrained on the class of permutation matrices, this 
framework enables us to find the optimal permutation matrix, 
i.e., subcarrier pairing, for relaying in both scenarios, without 
and with the direct path. The former is consistent with previous 
works, and the latter has never been shown before to the best of 
our knowledge. Furthermore, we prove that the optimal unitary 
linear processing matrix is in fact a permutation matrix which 
gives an optimal subcarrier pairing strategy. 

The rest of this paper is organized as follows. In Section HIl 
we present the system model and formulate the optimal linear 
processing matrix to maximize the achievable rate. Section Hill 
presents our analytical framework to determine the optimal 
permutation matrix and hence the optimal subcarrier pairing 
strategy. In Section HVl we show that subcarrier pairing in fact 
is the optimal linear processing strategy. We present simulation 
results to demonstrate the performance gain achieved through 
optimal subcarrier pairing in Section |V] and finally conclude 
in Section [Vl] 

II. Problem Statement 

A. System Model 

We consider a dual-hop relay network with a pair of source 
and destination nodes and a single AF relay node in an OFDM 
system with N subcarriers. We constrain ourselves to half- 
duplex transmission, where a relay node is either in trans- 
mission or reception but not simultaneously. The cooperative 
transmission takes place in two phases. In the first phase the 
source sends data through N subcarriers to the relay and 



destination simultaneously (if the direct path is available). The 
relay then performs linear processing of the received signals 
over N subcarriers, and forwards the amplified version of the 
processed signals to the destination. 

We denote the channel gain over subcarrier k from source to 
relay, from relay to destination, and from source to destination 
by /life, /i2fe, and /ipfc^ respectively, and Sfe the source signal 
transmitted on subcarrier k with power coefficient dsk- The 
received signals at the relay and destination in the first phase 
are given by 

y., = HiDsS + rir, (1) 
y^'^ = HoD,s + n(^\ (2) 

where = [yri, ■■■ , VthV and y^^^ = [y^^^ • • • , y^];^^ are 
the received signal vector at relay and destination, respectively. 
Hi = diag(/iii, • • • , hiN) and Ho = diag(/ioi, • • • , /iqat) are 
the corresponding channel matrices, and Ds = diag(ds) with 
ds = [dsi, • • • , dsNV being the power coefficient vector. The 

signals are i.i.d. with average unit power E[ss^] = I. More- 

r iT A (1) r (1) (i)iT 

over, rir = Ki, • • • , Tirjvj and = ["■di > ' ' ' ' "-dwJ are 

AWGN at the relay and the destination, with C7V(0, cr^I) 

and n|^^' ^ CN{0, tr^I), respectively. 

In the second phase, the received signal y^ is linearly 

processed at the relay with a unitary matrix W, and the relay 

retransmits the amplified version of the processed signals. The 

received signal vector at the destination is given by 

yf = H2D,Wy, + nf 

= H2D,W(HiD,s + n,) + 4'^ (3) 

where y^^^ = [y^^^ , ■ ■ ■ ,yf^Y , H2 = diag(/i2i, • • • , h2N f, 
and n^") = [n^dl ' ' ' ^ ^dN? ^ CN[^Ji,all). The power 
coefficient vector for the processed signal at the relay is 
denoted as d,., and we have D,. = diag(dr). Note that, W is 
a unitary matrix with WW-^ = I. 

Let Ps and Pr be the maximum average total power at 
the source and the relay, respectively. Thus, Ds must satisfy 
E||Dss|p = ||ds|p < Ps- Similarly, the linear processing 
matrix W and the power coefficient matrix at the relay 
must satisfy 

E||D,W(DsHis)||2+E||D,Wni||2 < P,, (4) 

To focus on the effect of the processing matrix W on the 
relay performance, in this study, we assume a pre-determined 
power allocation over subcarriers at the source, i.e., is 
given. At the relay, we assume that the relay equally amplifies 
the processed signal over the subcarriers, i.e., is a scalar 
of identity matrix = drl with 

/ 

dr = , rr (5) 

which is obtained based on Note that since W is a unitary 
matrix, it does not appeared in (|5]). 

In a relay system, the direct path between the source and 
destination may or maynot be available. When the direct path 



is available, we assume the receiver uses the maximum ratio 
combining to improve the reception and maximize the received 
SNR. 

B. Linear Processing and Achievable Rate 

We consider the achievable rate in such AF relay OFDM 
system. Regardless whether the direct path is available, we 
can rewrite the end-to-end system equation in the following 
general form 

y' = He,(W)s + ne, (6) 

where Heq(W) is the equivalent channel matrix and is a 
function of the processing matrix W, and the equivalent 
noise term. Then the system achievable rate is given by 

C(W) = i logdet(I + R-iHe,(W)Hf^(W)) (7) 

where R„ = E[negn^] is the covariance matrix of the 
equivalent noise term. The factor 1/2 reflects the half-duplex 
operation. 

For the conventional OFDM relaying without linear process- 
ing, i.e., W — I, the relay simply forwards the amplified signal 
to the destination over the same subcarrier However, such 
forwarding is in general not optimal in terms of maximizing 
the achievable rate. A special class of W is the permutation 
matrix H, for which linear processing reduces to subcarrier 
pairing. Such a scheme would uniquely couple a subcarrier 
over the first hop with a possibly different subcarrier over the 
second hop for signal relaying. This technique was studied 
recently in a few specific relay models y], Q, 16 1 and was 
shown to improve the rate. However, the optimality of such 
an approach remained unknown. 

Our goal in this study is to find the optimal linear processing 
matrix W* to maximize the achievable rate 

W*= argmax i logdet(I + R-iHe,(W)Hf,(W)). (8) 

In the following, we will first focus on the class of permutation 
matrices for relaying with and without direct path. We will 
then discuss the optimal W* for these scenarios. 

III. Linear Processing under Permutation: 
Subcarrier Pairing 

To solve (O, we first focus on the class of permutation 
matrices W H, and propose a general framework to find 
the optimal subcarrier pairing. This framework relies on the 
following result 

Lemma 1: Let P and Q be two diagonal matrices. For 

max det(I+ (PnQ)^(PHQ)), (9) 

among all possible permutation matrices H, the optimal H* 
is the one that maps the sorted absolute values of the diagonal 
entries of P to the sorted absolute value of the diagonal entries 
of Q. 

To show the above, we see that the objective function in ^ 
can essentially be rewritten as nr=i(-'^ + IPiPlftP)' where pi 



and qi are respectively the diagonal elements of P and Q. Let 
{\P(i) P} ^iid {\q[i) P} be the corresponding ordered sequences. 
Then it is not difficult to show that n"=i(l + IftPkiP) < 

nr=i(i + bwn9wH- 

A. Optimal Pairing for Relay without Direct Path 

In this case, the destination is out of the transmission zone 
of the source. From Q, the equivalent channel matrix, noise 
vector, and its covariance matrix in (|6]) in this case are given 
by 



He,(n) = H2D,nHiD, 
ne„ = H2D, Jlrir + 1 



H 



,(2) 



(10) 



Note that, except 11, all matrices in ( fTOb are diagonal. Using 
the property of the determinant, dct(I + AB) = det(I + BA), 
from (|7l), we can write the end-to-end achievable rate as 

c(n) = iiogdet(i + He,(n)^R;^iHeg(n)) (ii) 

Inserting the expression of R„, (fTTT l can be rewritten as 



c(n) = -x 



log det(I 



(R„ n HiD,)^ R„ nHiD,). 

P Q P Q 

(12) 



The i-th diagonal entries pi and qi of P and Q are pi = 
I 2 12 9' hiidsi, respectively. Following 

Lemma [1] the optimal 11 is to pair the ordered sequences 
of {\pi\} and or equivalently, to pairing the following 

two ordered sequences 



{SNR,,,,}, 



1 + a2 SNR^d, 



where SNR, 



_ \hi 



and SNRrd. = ^ 
the received SNR from source to relay, and from relay to 
destination, over the ith subcarrier, respectively. 



Note that, because f{x) 



is monotonically in- 



creasing for a > 0, ordering {\piW is equivalent to ordering 
\h2idri\. Under the noise-free relay assumption, i.e., a"^. = 0, 
the author of 13 proved that the sorted subcarrier pairing, 
which couples the relay subcarriers {|ciri/i2i|} with the trans- 
mit subcarriers {|c?si^ii|}, is the optimal pairing scheme. Our 
result shows that this pairing scheme is optimal for noisy 
relaying as well. 

B. Optimal Pairing Relay with Full Diversity 

We now consider the case when the direct path is available, 
i.e., relaying with full diversity. Through the maximum ratio 
combining technique, the observations from the direct path 
and the relay node can be coherently added. Define T2 = 



R„ "H2D^, Yi = HiD„ and To = ;f HqD,, where R„ is 
given in (UOl ). The achievable rate in this case is given by 

CiU) - i logdct (l + (T2nYi)^T2nYi + T^Yo 

'(13) 

To find the optimum 11 to maximize C(n), we again apply the 
result in Lemma [T] and all we need is to find the equivalent 
P and Q to express ( fT3] l as the form in (|9]). To do so, we 
re-arrange ( fT3] l as 

C(n) = ilogdet((I + T«Yo)- 



(I + + T^Yo)-i(T2nYi)^T2nTi 
= ilogdet(I + T?To) + 

i logdet (l + (T2nTi)(I + Y^To)-i(T2nTi)^ 



(14) 

where the second term of (fl4] l follows from the property 
det(I + AB) = dct(I + BA). Since the first term of (O 
is independent of IT, we are only interested in the second 
term as a function of 11, which can be written as 

C2(n) logdet (l + T2nTi(I + T^To)-^ 



(15) 



BA), we can set 
and ( fTSl l can then 



(T2nYi(I + T|^To) 

Again using dot (I + AB) = dot (I 
P = T2 and Q = Ti(I + T^Yo)- 
be transformed into the form of (|9]l. Based on this, we obtain 
the optimal subcarrier pairing scheme for relaying with full 
diversity. Examining the diagonal entries of P and Q, we 
conclude that the optimal pairing is essentially to order the 
following quantities over the input and output subcarriers of 
the relay, respectively. 



SNR,.,, 
1 + SNR.d, 



1 + ct2 SNR^d 



(16) 



are where SNR, 



, SNR^d^j 



and 

SNRsrf j — ' ° , are the received SNR from source to relay, 
from relay to destination, and from source to destination, over 
the ith subcarrier, respectively. Again, sorting { i+j^a^R.' ^ } 
is equivalent to sorting {SNRrd.i}- Thus, the optimal pairing 
for relaying with full diversity has a very clear strategy based 
on the SNR on each path: it is to match incoming and outgoing 
subcarriers of the relay, according to the SNR strength on 
the relay-destination subcarrier, and the relative ratio of SNR 
strengths on the source-relay and source-destination paths. 
Various sorting algorithms can be employed at the relay node 
with the computational complexity of 0{N \og2 N) [7 |. 

IV. The Optimal Linear Processing 

In previous section, we have considered the special class 
of permutation for linear processing, and have obtained the 
optimal permutation 11* in different relay scenarios using the 



proposed framework. Back to the general problem in the 
following result shows that the optimal permutation 11* is in 
fact the optimal linear processing matrix W*. 

Theorem 1: Let P and Q be two diagonal matrices. The 
solution to the following maximization 



max det (I - 

W:W«W=I 



(PWQ)^(PWQ)) (17) 



is W* = n*, where 11* is the solution in Lemma [T] 
Proof: See Appendix lAl 

Combining Theorem [T] and the results in Section |III1 the 
following result immediately follows: 

Corollary 1: The optimal W* for ^ in the relaying prob- 
lem with or without direct path considered in Section |III] is 

w* = n*. 

In other words, the optimal Unear processing is essentially 
subcarrier pairing with the optimal pairing strategy demon- 
strated in Section Hill 

V. Simulation Results 

In this section, we compare the performance of the opti- 
mal subcarrier pairing (SP) scheme with other non optimal 
schemes through Monte-Carlo simulations. We consider the 
scenario of relaying with direct path. We assume a 5MHz 
OFDM system with N = 128. The achievable rate is averaged 
over randomly generated multi-tap frequency selective chan- 
nels (with fixed channel order L = 11). A source-destination 
pair is placed at a distance dgd apart, and the distances between 
source and relay, and relay and destination are set at d^^ 
and drd, respectively. The pathloss exponent of 2 is assumed. 
We assume Pg — Pr, and Pg is equally allocated across 
subcarriers, i.e., d^j = 



We denote SNR = 

N Na 



as 



the average per subcarrier received SNR over the direct path. 

We first consider dsd = 20m, dsr — 6m, and drd = 16m. 
FigHldepicts the performance of the average rate per subcarrier 
vs. SNR, averaged over random channel realizations under the 
following four schemes: optimal SP scheme 11*; no SP used; a 
random W used; using IT* that is obtained assuming no direct 
path (i.e., (fT2li). The reason we consider the fourth scheme is 
that, in some cases, it may be easier for the relay to compute 
the optimal SP only based on the SNRs obtained on the two 
relay paths, although the receiver may use direct path signals 
for processing. We see that the optimal SP scheme outperforms 
all other schemes. When compared with the non SP scheme, 
it provides about IdB gain. We expect to see a substantial 
gap between the optimal SP and other schemes in the highly 
asymmetric and frequency-selective channels. 

Next, we study when subcarrier pairing is the most benefi- 
cial. In Fig|2] we show the effect of the relative SNR strengths, 
between relay and direct paths, on the achievable rate under 
different linear processing schemes. With fixed dgd = 20m, 
we vary the relay position between source and destination. 
We plot the average rate vs. the relative distance dsr /drd for 
SNR = 14dB. This figure shows that the performance of all 
schemes coincide, when relay is very close to destination, i.e.. 
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Fig. 1: Relaying with direct path: Rate vs. SNR 
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Fig. 2: Relaying with direct path: Rate vs. 
relative distance ^ 



SNRrd is very high. In this case, SP does not provide much 
performance gain. On the other hand, when the relay moves 
closer to source, the gain of using the optimal SP over other 
schemes become more substantial. 

VI. Conclusion 

In this paper, AF relaying with unitary linear processing at 
the relay in an OFDM system is considered. We proposed 
a general framework to analyze how to select the unitary 
processing matrix to maximize the system achievable rate. We 
first considered the class of permutation matrices, and based 
on the proposed framework, we derived the corresponding 
optimal permutation matrix, or subcarrier pairing scheme, for 
relaying scenarios with and without direct path for diversity. 
We further show the optimality of the so obtained permutation 
matrix among all unitary linear processing matrices for system 



achievable rate maximization. Simulation results also demon- 
strate the gain can be achieved through optimal subcarrier 
pairing as compared to the non-optimal hnear processing and 
non-pairing cases. 

Appendix A 
Proof of Theorem[T| 

We prove this by induction. We provide a brief description 
of the steps leading to the conclusion. 

For = 2, we need to show that the permutation matrix 
obtained from sorted subcarrier pairing is optimal. In order 
to do so, we first parameterize a 2 x 2 unitary matrix and 
then find its optimal entries. From the general parametrization 
of unitary matrices flS), the 2x2 unitary matrix W can be 
represented by 4 independent parameters 



W 











cosf 
sin^ 



— smf 
cos 9 





aj4>3 



(18) 



Substituting ( fTsT i in (fTTT i and taking derivative of the determi- 
nant with respect to four existing variables and 6 of 
W, we reahze that the optimal 6* = 0, when pi > p2 and 
qi > q2, or pi > p2 and qi > q2- The other variables are 
canceled out in the determinant. This observation shows that 
in the case of iV = 2, W* is a permutation matrix. 

Assume this is true for N = n — 1. For N = n. Let 
A„ = (P„W„Q„)^(P„W„Q„), where subscript n denotes 
the matrix dimension. Since (I„ + A„) is a positive definite 
matrix, it has the following property f9\ 

det(I„ + A„) < (1 + a„„) dct(I„_i + A„_i) (19) 

where a„„ ~ (A„)„„. Re-arranging the expression of the 
determinant, we have 

det(I„ + A„) = det(I„ + W^Q^QWPP^) (20) 

Let W = (wi ••• w„), P = diag(p) with p = 
[pi,-- ■ ,PnV, and Q = diag((7i, • • • , q„). We have 



-^71 



|p|'0Wi)^|gi|2wi 



IpP w„)^q?wi 



(p2 0wi)^|(7„pw„ 



|pp0w„)^|(7„|2w„ 



A„_i IpP [wi • • • w„„i]^|(j„pw„ 

PP w,f |qi|2W„_i (|p|2 w„)^|(7„|2w„ 

(21) 

From ([T9j, we have 

71 

det(I„+A„) < {l+J2\P^\^\Win\^\qn\^)det{In-l+An-l) 

(22) 

Note that 



for any permutation matrix 11. Let \q(^n) P — niax{|gip}; then 
w.l.o.g., we let |q„p |g(ri)P- Thus 

n 

det(I„ + A„) < {l+^p^^\w^n\^ql)det{In-i+An-i) 

i 

< (1 + |p(„)ng(„)ndet(I„_i + A„_i) 

(24) 

where |p(„)|2 = max{|p,|2}. 

Since W*_]^ = H* _]^, by the Hardmard inequality, we have 



maxdet(I„_i+A„_i(W)) = [] {l + \p(^,)\%,^f)) (25) 

i 

where both p^^^ and g^^-j are sorted in ascending order Thus 

n-l 

det(i„ + A„) < {i + \Pin)\'hn)\') n (i + lml'kwn) 

i 

with equality if and only if w„ = [0, • • • , 0, 1, 0, • • • , 0]^ has 
entry 1 at the (ri)th position. This implies that 



(26) 



where 'W'^_^ is W„_i with an additional row of zeros 
inserted at the {n)th row. Hence, W„ is a permutation matrix. 
Therefore |I„ + A„|, i.e., the objective in (fTTI l is maximized 

by w; = n;,. ■ 
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